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Abstract 

We discuss the impact of the soft degrees of freedom inside the B c meson on its 
rate in the semi-leptonic decay B c — > X£pg where X denotes light hadrons below 
the D° threshold. In particular we identify contributions involving soft hadrons 
which are non- vanishing in the limit of massless leptons. These contributions become 
relevant for a measurement of the purely leptonic B c decay rate, which due to helicity 
suppression involves a factor m| and thus is much smaller than the contributions 
involving soft hadrons. 
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1 Introduction 



Within the standard model of elementary particle physics the B c is the only meson which 
contains a heavy quark- ant iquark (QQ) pair of different flavour. After more then two 
decades of purely theoretical investigations its observation by the CDF collaboration in 
1998 pushed open the gate for experimental exploration. The newest and forthcoming 
generation of experiments (Tevatron and LHC) will be able to determine the B c properties 
that are hardly known at present 0. 

Theoretically, the fact that the B c is a bound state of heavy quarks offers an interesting 
possibility for an efficient way of calculating its properties. While the bound state is 
dominated by soft scales the heavy quark mass m Q which is the typical scale for production 
and decay processes of such mesons is much larger. Due to this hierarchy of scales it is 
obvious to use an effective theory approach. One could treat the B c within the Heavy 
Quark Effective Theory (HQET) at intermediate scales (i.e. where the charm (anti)quark 
is considered to be light) but this does not take into account that m c can still be considered 
a perturbative scale. Assuming that the c quark is also heavy (i.e. using that both m b , m c 3> 
Aqcd) we may work with approaches developed for quarkonia-like systems. In this case the 
B c stands somehow in between of charmonia and bottomonia which usually are investigated 
within the framework of Non-relativistic Quantumchromodynamics (NRQCD) Q. 

In this article we concentrate on the semi-leptonic B c decay B c — > Xiv^ where X denotes 
only light hadrons which means hadrons below the charm threshold. In this kinematic 
region the B c has to decay through electroweak annihilation. The light hadrons thus 
originate from gluons or other light degrees of freedom inside the B c meson which makes 
this decay mode particularly interesting, since it tests the light degrees of freedom in a 
quarkonia-like system. But also for practical purposes it is of relevance: We will show that 
the semi-leptonic decay rate exceeds the leptonic one (at least for light leptons) even in the 
limit of extremely soft light hadrons which could complicate substantially the measurement 
of the latter. 

The total rate of the purely leptonic mode is given by the expression 
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where Gf is the Fermi constant, V c b the CKM matrix element corresponding to the un- 
derlying partonic decay process be — ► tv^ and M is the mass of the B c . The B c decay 
constant f& is defined via 



if B y:= (0\h^c\B c (p)) 



(2) 



Obviously the total rate (ffl) is suppressed by the ratio mj/M 2 of the lepton and the B c 
mass squared. This is caused by the helicity flip which is necessary in a back-to-back decay 
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of a spinless particle in the ground state into a left-handed lepton (right-handed antilepton) 
and its right-handed antineutrino (left-handed neutrino). 

Contrary to the quark model where the (be) systems has to be in a 0~ state, NRQCD and 
related effective approaches allow for higher Fock states, i.e. the (be) pair also can exist with 
other quantum numbers than J p = 0~ inside the B c . Although these states are suppressed 
by non-perturbative factors as v 2 in NRQCD[] their contribution becomes measurable if the 
mj suppression of the leading Fock state is stronger than the v 2 suppression of a subleading 
Fock state which does not suffer from helicity suppression. However, these contributions 
are not purely leptonic anymore since in such a case the soft degrees of freedom inside the 
higher Fock states should give rise to some soft hadronic signal in the final state. 

Thus we expect, for instance, semi-leptonic contributions related to the NRQCD matrix 
elements (Oi( 3 Si)) and (Oi^Pi)) where the short distance coefficients arise without any 
m 2 /M 2 suppression since it is not necessary to flip any helicity if the decaying (be) pair 
has spin one or is in a P wave state. Additionally there also should be unsuppressed 
contributions from higher Fock states with the (be) pair in a 1 Sq state since the soft degrees 
of freedom inside the B c can carry away some angular momentum or spin to prevent 
the helicity suppression in the leptonic sector. Thus the corresponding short distance 
coefficients should (at least for light leptons) overcome the enhancement of the leading 
order matrix element (Oi^Sq)) against the higher orders in the non-relativistic expansion. 
They also should exceed the short distance coefficients from radiative leptonic B c decays 
where the prize of avoiding helicity suppression by radiating off an additional photon is 
paid for by a factor a/ (47r). 

In order to calculate the lepton energy spectrum in the decay B c — > Xlv^ we apply an 
effective approach almost equivalent to standard NRQCD. Even though we use a different 
notation namely a covariant one borrowed from HQET our effective Lagrangian matches 
the NRQCD one. However, we enlarge the operator basis at dimension-8 by operators 
related with the centre-of-mass (cms) momentum of the QQ pair inside the bound state 
since we expect contributions that are not helicity suppressed just from these operators. 
According to the standard NRQCD velocity rules the matrix elements of such operators are 
suppressed with respect to matrix elements of dimension-8 operators that merely depend 
on the relative momentum of quark and antiquark. On the other hand, it is known that 
these power counting rules are only applicable for systems which obey the scale hierarchy 
m Q > m Q v > m Q v 2 ~ A QCD 0. 

Fleming et al. |J proposed to keep the standard NRQCD velocity rules only for bottomonia 
while they suggest to switch to a dimensional power counting in the charmonium sector. 
The gist of this analysis is a new weighting of a spin-flip transition that is compatible 
with data M and possibly could explain some of the NRQCD drawbacks concerning the 

■""We have chosen v instead of the commonly used v to denote the typical (average) non-relativistic 
(anti)quark velocity inside the quarkonia, i.e. more precisely the absolute value of its three vector, to avoid 
confusion with the covariant four vector we need to define the quarkonia velocity. 
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J /if) polarization at high transverse momentum and the endpoint spectrum in radiative 
quarkonia decays. 

Here we do not only discuss the relative weight of a dipole and a spin-flip transition. 
Rather the semi-leptonic decay rate may also help to clarify the power counting of the 
cms derivatives on the QQ pair. Our procedure follows the strategy of Heavy Quarkonia 
Effective Theory (HQQET) P,^. We just eliminate some redundancies and rotate the 
operator basis into the familiar one of NRQCD. 

In the following we will start with setting up the basics of the effective approach. We sketch 
how the Lagrangian is constructed and how B c annihilation processes are implemented. 
Furthermore we provide the basis of four fermion operators up to dimension-8 and discuss 
the power counting rules. After applying this framework on the lepton energy spectrum 
in the decay B c — > Xti>i we address the structure of its endpoint region. We also relate 
our result to the shape function improved one of NRQCD. Finally we discuss the total 
semi-leptonic width of the B c and implied consequences on measuring its purely leptonic 
width. 



2 Effective field theory for the B c 

In this section we shall briefly discuss the effective field theory set up as it is appropriate for 
the case of a B c decay. In general, quarkonia are more complicated than mesons containing 
only one heavy (anti) quark due to the fact that several scales are involved, some of which 
are dynamically generated. 



2.1 Lagrangian and fields 

In the first step one uses the fact that the heavy quark mass is large compared to any other 
scale in the problem. Thus one performs a 1/m.Q expansion of the full QCD Lagrangian, 
which can be done either by integrating out the heavy degrees of freedom, which are the 
small components of the heavy quark spinor |I| , or by performing a sequence of Foldy- 
Wouthuysen transformations ||11|| . In both cases one obtains a l/rriq expansion of the 



heavy quark field as well as of the Lagrangian. Although these expansions look slightly 
different, the results for the matrix elements are the same. Up to I/^q one obtains 

C v = fc« (ivD)h^ + ^—{K v + M v ) + O ( \ ) (3) 

Zm Q \ m Q J 

with 

K v = h^iD^g^ - Vfl v v )iD v h^ , (4a) 
M v = !|i/4+)(-^)G M ^« (4b) 
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where = P^/M is the four velocity of the heavy meson and h v the field operator of a 
heavy quark in the limit m Q — > oo. The terms K v and M v correspond to the kinetic energy 
and chromo magnetic moment, respectively, where the gluonic field strength tensor G^ v is 
defined by the relation 

igsG^u := [iD^, iD u ] . (5) 



Likewise, we obtain the corresponding relations for an antiquark in the infinite mass limit 
by reversing the sign of the velocity v — > — v. 

In order to describe a quarkonium-like system one wants to treat both the quark and the 
antiquark in the infinite mass limit, assuming that both are sitting in a bound state moving 
with a velocity v and that both heavy constituents have small residual momentum. Thus 
one would start with the Lagrangian 

£ s tat. = hi + \tvD)h^ _ h$\iv'D)htf (6) 

and add the corresponding pieces of the Lagrangian for the light degrees of freedom. 

However, it is known that this simple ansatz misses an important piece of physics, which 
is that the dynamics responsible for the binding generates additional scales. These scales 
are the size and the binding energy of the bound system and are related to the relative 
(three-)velocity t> r ei which can be written as 

^rel= vV') 2 -l- (?) 



Clearly in a quarkonium we have t> re i -C 1 and in the limit v — > v' the Lagrangian ([]) 
exhibits certain pathologies indicating that bound states cannot be described in the static 
limit. As a example one may consider the matrix element (H stSit ] hv^h^ |0) in the static 
limit; computation of the anomalous dimension reveals an imaginary part in the anomalous 
dimension of the respective current Th[ 
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as v — ► v 

v' — > v leads to a phase factor in the solution of the renormalization group equation which is 
related to the coulombic part of the one-gluon exchange [|I3| , [!Hf 
of heavy bound states. 



The divergence in the limit 
ion group equation which is 
This reflects the possibility 



This problem can be fixed by including the kinetic term into the leading order Lagrangian, 
i.e. one switches to a non-relativistic description. The Lagrangian then reads 



C = ivD 



2m^ 



- Ur ] ivD + 



(iD 



2m^ 



where the component iD^ orthogonal the t> M is given by 

iD\ := (g^ - v"v v )iD v 



(9) 
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The form of the leading order term evidently shows that is impossible to operate with 
a purely dimensional power counting in this case, since the equation of motion is not 
homogeneous in l/m Q : 

(ivD)hl* = ±^htK (10) 

Usually it is argued that this corresponds to an expansion in the relative velocity v which 
is assumed to be small. In fact, the non-relativistic approximation (J8|) can generate bound 
states of (inverse) size m Q v and of binding energy m Q v 2 which introduce two more, dynam- 
ically generated scales into the problem. This fact makes the construction of an effective 
field theory approach more complicated. 

We also note that the chromomagnetic moment M v is taken to be higher order compared 
to the kinetic energy K v although it is of the same mass dimension. This means that 
the antisymmetric combination of two covariant derivatives acting on the field hv is 
suppressed with respect to the symmetric combination: 

[iD^iD^h^ < {iD^iD v }h^ . (11) 



The standard rules for power counting in non-relativistic system can be found in [15], but 
we shall return to this question at the end of this section. 

The static as well as the non-relativistic limit of QCD has (compared to full QCD) an 
additional spin symmetry which is evident from the Lagrangian (^|). In the case of a 
quarkonium-like system this symmetry predicts quadruplets which are degenerate up to 
terms of order v 2 . For the ground state quarkonia this quadruplet consists of the 0~ states 
and the three polarizations of the 1~ state. For excited quarkonia this quadruplet structure 
is more complicated and involves a possible orbital angular momentum. 

Finally one may formulate NRQCD in the rest frame of the quarkonium in which v = (1,0) 
which transforms @ into the usual NRQCD Lagrangian. Rewriting the (relativistic) four- 
dimensional fields hv in terms of the (non-relativistic) two-dimensional spinors ip and \ 
for quark and antiquark, respectively, 

one actually gets the NRQCD Lagrangian 

^4 A+ ^ + 4 A - 1 ?)*- (i3) 
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2.2 Operator Product Expansion for semileptonic B c decays 



The main contribution to semileptonic B c decays comes from the part of the effective 
Hamiltonian mediating b — > c transitions 

H cS = ^V cb (b L ^c L )(v e:L rh) (14) 

where the subscript L denotes the left-handed component of the the corresponding fermion. 
In order to consider the inclusive case we write the rate as 

dT = sc^jVctl y_A_ j d , xe _ iqx {Bcl ( - 6L(;r)7MCL(a;))(c - L( o) 7 ^ i( o)) \ Bc ) I^(q) (15) 

where I^^q) denotes the leptonic part. 

We perform an operator product expansion (OPE) for the product of the hadronic currents 
making use of the fact that both m c and m b are large scales, which can be removed from 
the matrix element by redefining the phase of the heavy fields. We consider the decay of 
a B~ which consists of a anti-c and a b quark; therefore we define 

b(x) =: e- im » vx B v {x) and c(x) =: e im * vx C v (x) (16) 

where B v and C v eventually become the nonrelativistic fields considered in the last para- 
graph, and v is the four velocity of the B c meson. In terms of these fields we have for the 
hadronic tensor 

= Jd 4 x exp{-t(q - M h - C v)x} (B c (v)\ (^ iL (x) 7/1 C„, L (o:))(a,L(0)7, J B„ )L (0)) \B c (v)) 

(17) 

where 

M b£ = m b + m c (18) 

is the sum of the quark masses. 

We perform an OPE assuming that (q — M b5 v) 2 and qv — M b5 are of order M b 2 g and M b5 , 
respectively. The differential rate takes the form 

dT = J2Uri(0[n]} (19) 

where C n are short distance coefficients which are calculated in perturbation theory and 
(0[n\) are forward matrix elements of local operators 0[n] with non-relativistic B c states. 
The parameter fi is the renormalization scale; the dependence of the coefficients on fi is 
canceled by the corresponding dependence of the matrix elements, such that the total rate 
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is independent of fi. The sum runs over an infinite set of operators which are of increasing 
mass dimension; this dimension is compensated by the coefficients C n which contain only 
the perturbative and large mass scale M& g , while the matrix elements (0[n]) do not depend 
on this mass scales any more; they depend only on a non-perturbative small mass scale A 
corresponding e.g. to the binding energy of a coulombic system. Thus higher dimensional 



operators are expected to be suppressed by powers of A/M b5 and thus the series in (|19]) 
may be truncated and we shall consider only the first nontrivial terms. 

In leading order we have two operators of dimension-6: 

0(ixi) t 1 ^] := - (h^^Ch^) (^H 75 C/i«) , (20a) 
0(ixi) flsf>] := - (h^^Ch^) {hl-hlCh^) . (20b) 

Here we have introduced the spectral notation n = 2S + 1 Ly^ where the upper index C — 1,8 
denote the the colour states that are indicated by C = 1 and C = T a for singlet and octet 
states, respectively, on the right hand side of equations (|20|). Using flT2| ) we recognize that 
the leading order four fermion operators again coincide with the one of NRQCD 

( ixi) P-Sf } ] = O c {%) , ( ixi) [ 3 S[ C) ] = O c (%) • (21) 

In subleading order our set of operators is more general than the standard NRQCD basis, 
since we also include contributions related to the motion of the QQ pair, more precisely its 
centre of mass, inside the quarkonium. It is known that these contributions are important 
in certain kinematical regions, e.g. in the endpoint of photon spectra in radiative quarkonia 
decays (|I^,|I7||, even though they are missed by naive use of NRQCD. To disentangle the 



relative motion of the quark and the antiquark from the motion of its centre of mass we 
decompose their residual momentum inside the meson according to 

• Residual cms momentum (RCM): 

iD, (Q(±> TC Qj*>) := g(±) TCiiD,) + gW (iD^TC Q?) (22) 

• Residual relative momentum (RRM): 

Qt ] r (0,) c gw := \ [g?) rc(iD,) gf) - g(±) (wjrc gf>] (23) 

Besides we introduce T^") := (T^ + T" At )/2 — g^T^/A to denote the symmetric traceless 
part of a tensor. In this way we can distinguish three types of dimension-8 operators: 
RCM x RCM: 

0(0x0)^°°] := " K +) 75^/i«)] [iD* {h^^Ch^)] (24a) 
0(CxC)[ 3 ^ C) ] := - [iDi {WriCht))] [iD» ± (hl-hlChM)] (24b) 
PfCxOpf'] := - [iDfr (h^^Ch^)] [iD^ (h^llCh^)] (24c) 
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RCM x RRM: 



C(cxr)[ 1 'S'o C ' ] 

^CxR)^] 



RRM x RRM: 

0(RxR)N C) ] : = 



iDi (hi +) 1*C [^") 75 (-l^j G 

«tf K +) 7^7 /£">)] [&H 7l (-IB^) C 

iic W)] W ii (-{d$i c W 



+ h.c. 
+ h.c. 
f h.c. 



75 



^ 75 



1 



(RxRJL^O 



4 TA-L 



Ch^ 



(25a) 
(25b) 
(25c) 

(26a) 
(26b) 



(26c) 
(26d) 



C(RxR)[ 3 -^2 ] 



7 ? (RxR)[5 , d 



^RxR^P] 
^RxR)^, 3 ^] 



^ +) 75C^ 



h^jsi-lD 1 ) Ch^ 



+ h.c 
+ h.c 



(27a) 
(27b) 



h.c. 



(27c) 



The main difference to the operator basis of NRQCD is the inclusion of operators of the 
type RCM x RCM and RCM x RRM. Dealing with the cms momentum of the QQ pair 
inside the quarkonium they were neglected in the original work [|J since they do not yield 
contributions from the dominant Fock state \QQ). Nevertheless it was shown that they 



play an important role in quarkonia decay JTj] and production |18[ processes as long as 
the observable is not inclusive enough. More precisely these operators are responsible 
for the shift from the partonic to the hadronic endpoint in energy spectra caused by the 
mass difference M — 2m Q ~ m Q v 2 . After resumming the leading cms operators to so-called 
shape functions the NRQCD prediction is extendible up to higher energy values. Numerical 
analyses of this effect using a model consistent with the shape function formalism indicate 
on the one hand a natural solution of the HERA problem in the J/ip photoproduction 



channel [|19],|2(| but generate on the other hand some new obscurities in the radiative T 
decay 0. 

The operators or the RRM x RRM sector common with NRQCD are related to its 
dimension-8 operators in the following way (cf. eq. ([12])): 



O 



(RxR)l 



|2S+1 L <,C)] 



a 



V {KxR) [ 2S+1 Lf] = Vc 



: 2S+1 lj) 



(2* 
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2.3 Power counting 



Instead of applying strictly the standard NRQCD velocity scaling rules we propose to use 
a slightly generalized power counting scheme. In particular we want to direct the readers 
attention on contributions related to the cms movement of the QQ pair inside the meson. 
These contributions potentially exceed the prediction by standard NRQCD power counting, 
which we will critically recapitulate first. 

To determine the suppression factor of NRQCD matrix elements of four-fermion operators 
we have to consider both the scaling of the operator itself and additional orders of v coming 
from the Fock decomposition of the hadronic state: 

\H) = E WM \QQ[n]) + cf [n] \QQ[n}g) + c%[n] \QQ[n]gg) + ...). (29) 

n 

While the standard velocity rules for the operators are straightforward the assignment of 
suppression factors to the coefficients Cj [n] in (^) is quite complicated since it depends on 
the non-perturbative dynamics inside the bound state. In the original NRQCD paper || 
quarkonia states are assumed to obey the scale hierarchy m Q ^> m Q v ^> m Q v 2 ~ Aqcd- In 
other words their binding is assumed to be coulombic which enables the use of NRQED 



velocity rules of fla| . In particular this implies the validity of multipole expansion for 
calculating the coefficients, i.e. a spin-flip transition (Ml) is suppressed by v 2 in contrast 
to a dipole transition (El) which only accounts for v. 

Beneke mentioned that this is no longer true if the scale hierarchy is m Q 3> m Q v ~ Aqcd 3> 
m Q v 2 ||. He proposed to consider Ml transitions as Xv where A = v for coulombic systems 
while A ~ 0(1) for m Q v ~ Aqcd- A recent analysis || elaborates this feature to claim that 
the standard NRQCD velocity rules are only applicable in the bottomonium sector while 
it suggests a HQET-like power counting scheme in Aqcd/^q for charmonia. 

We take a further step towards generalizing the counting rules questioning the scaling of a 
cms derivative inside a NRQCD operator. According to 

| Pe5 | = | Pjoft | = Em = ^1 = ^2(M - 2,n Q) (30) 

the total momentum Pqq of the QQ pair inside a quarkonium is connected with its binding 
energy M — 2mq. Therefore in the standard formalism the cms derivative on a quark- 
antiquark bilinear accounts as v 2 . However, if the binding of the quarkonia-like system is 
not coulombic anymore the binding energy no longer should be purely coulombic either. 
Rather there also should be non- negligible contributions of order Aqcd/^q ^ v 2 . Hence 
we assume a cms derivative to be of order kv keeping in mind that the standard NRQCD 
power counting is restored by k = v. 

With this set of rules the B c decay matrix elements of the operators up to dimension-8 
scale as indicated in Fig. |T]. Note that colour selection rules forbid the existence of solely 
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one gluon in the higher Fock state if we restrict the perturbative part to leading order a s . 
To estimate the leading contributions according to the non-relativistic expansion we have 
to know if the quarkonia-like B c behaves rather as bottomonia or as charmonia. In the 
first case (k, A — > v) the first five configurations of Fig. [I] would contribute up to order 
0(v 8 ) while for k, A — > 1 the main matrix elements are the ones of the first row and the 
first column in Fig. [|. 



3 Lepton energy spectrum 

Using these assumptions we shall now calculate the lepton energy spectrum dT/dz in the 
semi-leptonic decay B c — > Xlj?i where X denote only light hadrons. 

z ' = w K < 31) 

is the lepton energy Eg in the B c rest frame normalized on half of the sum of the heavy 
quark masses. 

Note that Mbc is slightly smaller than the mass M of the B c . The mass difference M — 
just count for the binding energy of the (be) system inside the B c meson which is of the 
order v 2 within usual NRQCD. However, this has to be taken into account to get the 



correct endpoint of the energy spectrum \T?\. We will return to this point later again. 



The underlying partonic process is the electroweak decay b + c — > tvg. Thus we get 

<l[ ' ' ( ] 2 \V cb \ 2 [d'x £L e -M«-M b M Uq) {Bcl j, {x) m w . (32) 



dz M \ V2 J J (2tt) 

The B c rest frame is characterized by — (1, 0). There the tensor solely depends on 
the momentum q that flow through the effective vertex 

dp e dp De (2n) A 5 A (q -p e - pp e ) 8{z - — ) (33) 

which means simultaneously that only the symmetric terms of the leptonic tensor 

IT = 2 (rffi t + pU t + <T \Vi ■ P* t ) - ie^plvQ (34) 

survive. Let us turn to the hadronic matrix element now. Since we have factorized off 
their main phase space dependence the fields B v (x) and C v (x) in the left-handed current 

J${x) = B v {x) 1>1{l ~ lb) C v {x) (35) 
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only depend on residual momenta. Hence we can expand the hadronic matrix element in 
the coordinate space around x = after we have ensured hermiticity by symmetrizing the 
matrix element previously: 

(B c \ JZ{x)%(0) \B C ) -> (B c \ i[J£(s)#(0) + JZ(0)r v (-x)} \B C ) =: . (36) 

Then the expansion in x correspond to an expansion in the residual cms momentum of the 
(fee) pair inside the B c . Up to second order we get 

JC(*) = iC(0) + iS-iCC^U ("««) + 5 i^i^iC^U H^)H^) + • • • (37) 

where the currents on their parts are also expansions parameterizing the residual relative 
momentum of the (fee) system. It is given in eq. (|66"D of the Appendix. 

The x dependence in the phase space integrals of the form 

d 4 x -0^ {1, (-ix a ), {-ix a ){-ix p )} e -^- M »^ V(g) (38) 

is regarded by expressing the factor — ix a as derivative d/dq a on the exponential and 
subsequent integration by parts. At this stage arise higher derivative of delta functions, 
and the result for the spectrum is of the structure 

U = T ^AS{1 + e-z) + — [B 6{1 + e-z)+ B 1 5'(l + s - z) + B 2 6"(l +£-z)}\. 

(39) 

The fact that the differential decay rate is proportional to localized distributions (i.e. S 
functions and its derivatives) reflects the two body structure of the decay at (partonic) 
tree level; the kinematics fix the lepton energy to be z = 1 + e where 



(40) 



is the rescaled lepton mass squared. Clearly e is a small quantity for the electron and the 
muon case, in which we can use only the leading term. However, for the r the parameter 
e is sizable and we thus keep the full dependence on e. 

The coefficient A contains contributions from both the spin singlet and the spin triplet 
dimension-6 operator: 

A = Aixi) t 1 ^] (B c \ ( ixi) t 1 ^] \B C ) + A {lxl) M^] (B c \ Cdxi)^] \B e ) . (41) 

We expect that the leading term represented by ^4.(ixi)[ 1 'S'o ] ^ s suffering from helicity 
suppression just like the purely leptonic decay. In the inclusive case this is due to spin 
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symmetry making its contribution to R^ v entirely proportional to v^v u . If vacuum insertion 
is applied, this contribution becomes exactly equal to the one from the purely leptonic case. 
In fact r ^4(ixi) PSo^] reproduces the total rate ([1]) of the purely leptonic decay: Defining 

r - Gl l V cb \ 2 M? s (42) 



8ttM 

the short distance coefficient reads 

Aixi)^] = e(l - e) 2 . (43) 

While the Wilson coefficient of the spin singlet operator explicitly shows a factor e we do 
not expect any helicity suppression in the spin triplet case: We obtain for the corresponding 
short distance coefficient 

AiM)m = (i^m±£) (44, 

which does not vanish as e — > 0. From this observation we deduce that the spin triplet 
operator could become important even though its NRQCD matrix element is suppressed 
by X 2 v 6 with respect to the leading matrix element of the operator O/xxi^Sq ]. If we 
compare the magnitude of the suppression factors 

AixiifgQ ^ A 2- 6 „ m o (lxl) [*s?]\b c ) 

AixDt 3 ^] ~ 2 ~ ~ (B c \ Odxi,^] |S C ) 

we realize that the spin triplet contribution overcomes the spin singlet contribution if the 
lepton mass is small enough. For sure this condition is fulfilled in the electron case. For 
muons both contributions are comparable although it depends on the value of A and v for 
B c . In the case £ = r we can neglect the spin triplet part. 

The coefficients Bi arise with a factor l/(2m) 2 where m is the effective (= twice the 
reduced) mass of the b quark and the c antiquark. Furthermore we have introduced the 
dimensionless variable 77 to parameterize their mass difference: 

2m b m c m b -m c 

m := , 77 := . (46) 

rrib + m c m b + m c 

Generally operators of higher powers in the inverse mass m are accounted to be of higher 
order. However, there is a big caveat in the endpoint region of the energy spectrum. 
Because the expansion is an expansion in terms of 1/ (Mf, 5 — 2Eg) rather than in terms of 
l/(2m) where 

2m = (1 - V 2 )M bs (47) 

it breaks down in the endpoint region. In the framework of NRQCD it was shown that 
the availability of the theoretical prediction can be extended to higher energy values if the 
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contributions of the order [m Q v 2 /(m Q — Ee)] n are summed up to so-called shape functions 
T7j. Thereby the endpoint also shifts from the partonic value characterized by mass of 



the QQ pair to the hadronic value which is a function of the mass of the decaying particle. 
In our case higher order derivatives of the leading delta function 5(1 + e — z) represent 
contributions to this shape function. We will elaborate this in Subsection |3.1| more detailed. 



First we present the result for the coefficients Bi in (|39|). Due to heavy quark spin symmetry 
the spin triplet operators get an additional suppression factor compared to the spin singlet 
operators when they are inserted into the B c matrix element. On the other hand they 
should not suffer from helicity suppression as we have seen from On x i) 

[ 3 S[ 1] ]. Therefore 

we will not only keep the spin singlet contributions to Bq but also the most relevant triplet 
operator 0(r x r) [ 3 Pj] according to Fig. |l|: 



°0 - C ( C xC) 




(fiel^Cxot^] 


\B C }~ 


, «(0) 
r °(CxR) 


['si 11 


(B c \ C(CxR) 




\B C ) 


+ C (0) f 

^ L "(RxR)i 




(ficl^RxR)^] 




- B {0) \ 




(B c \ C(RxR) 




\B C ) 


^ °(RxR) 


[3P« 


(S c |0(R X R)[ 3 Po (1) 


\B e ) 













(48) 



For light leptons the pieces which do no suffer from helicity suppression can be as sizable as 
the dimension-6 contribution A. If we consider the limit e — > approximately fulfilled by 
electrons we should find non- vanishing contributions from ^(r x rJ 1 Pi ] and $(r x r)[ 3 -Pq ] 
where the quark-antiquark pair is in a P-wave state. Additionally we expect the coefficients 
^(CxClt^o ] an d ^[cxRjt 1 ^ 1 '] to survive in this limit. There the soft degrees of freedom 
moving relative to the quark-antiquark pair carry away some spin or angular momentum. 
Indeed WrxrJ^o ] i s the om y pi ece of Bq that is proportional to e: 



»(0) rlo(l)l 
°(CxC)l°0 J 



»(0) rlo(l)l 
°(CxR)L°0 J 

U (RxR)l°0 J 

B^l^pV] 



2(6 - 8r] 2 + 3r] 4 



8r](2 - 



8T]£ 



I + r] 2 + 2(3 - 6r] 2 + 2r] 4 )e 



1 - e H -e 2 



(RxR)l 

o(0) r3p(l)i 
D (RxR)l r 



-As [l - 2(2 - 7] 2 )e + (3 - 2r] 2 )e 2 ] , 
8 - As [3(1 - rj 2 ) + 6t] 2 e - (1 + 3rj 2 )e 2 ] . 



(49a) 

(49b) 
(49c) 
(49d) 
(49e) 



Note that the terms odd in 77 changes their sign if we consider P+ — > X£ + i>£ instead of 
B~ — > X£i?£. However, this does not mean that the result violates CP invariance, since the 
sign of the operator O^cxR)^^} also changes. Furthermore we have used heavy quark 
spin symmetry to summarize the P wave channels 3 Pq, 3 P\, and 3 P2- 



(B c \ C(RxR)[ 3 Pj 1} ] \B C ) = (2J+ 1) (B c \ (RxR) [ 3 P o (1) ] \B C ) 



(50) 
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Next we deal with the contributions B\ and £> 2 . As mentioned above they become relevant 
close to the endpoint since they are more singular than the leading contribution. Once the 
S function and its derivatives are resummed we will find a shape function enhancement 
close to the endpoint. For reasons which will become clear in a moment we also keep the 
£>i pieces stemming from C\cxC) [ 3 S[^] and V(cxR)[ 3 S[^]: 

B, = B^S^] (B c \ ( cxo I 1 ^] W + BgxR)^ 1 '] (B c \ 0(cxR)[ 1 sf ) ] \B e ) 

+ ^R)! 1 *'] W ^xRlt 1 ^] \B C ) + Bg xC) pSf >] (S c | ^cxof 3 ^] |5 C > (51) 
+ BgxR)^] < 5 d ©(cxr)^] \B C ) + Cg> xR) pSf >] (5 C | Pcrxr)^ 1 )] \B C ) . 

While the spin singlet contributions are still helicity suppressed the coefficients correspond- 
ing to a quark-antiquark pair in a spin triplet state do not vanish for e — > 0: 

Sfcxo^] = ^ £ [ 13 " " 3 ( 7 " 6 ^ + 3^ 2 - (11 - IOtjV] , (52a) 



B^p^ = - 2 -^ Y ^e[l-e]\ (52b) 



CgxR)! 1 ^] = "2(1- [l-e] S , (52c) 

Bfgxo^l = + ^[29 + V 2 - 15(1 + ^ (52d) 

- (17 - 27r] 2 )e 2 - 15(1 - r] 2 )e 3 ] , 

= - M f^ + [" " * + + , (IB.) 

C,^"] = + [5 - 3, - ^ + , 3 ] ■ (52f) 

Finally we quote the results for the coefficient B^ 2 \ In that case we only find contributions 
to ROM x RCM operators: 

B 2 = B^SU] (B c \ ^cxotf] \B C ) + B$ xC) flSf >] (S c | ^cxC)^] |S C > 

+ Cg ) xC) ps{ 1) ] (Scl^cxcjpS^llflc) 
where 

Bg^^fi-^sll-^, (54a) 
BfSxC, P^l = Ii ^7r l! " [25 - 44. + 38. 2 - 9, 3 - 3. 4 ] , (54b) 

cfio = + [ s - ^ + ^ - ^ + £ i • ( 54c ) 

Again only the spin singlet pieces show helicity suppression. However, this feature of the 
coefficients of higher order derivatives of the delta function does not surprise since it is given 
by the phase space structure of the process as we will explain in the following subsection. 



(53) 
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3.1 Shape function formalism 



The spectrum within the operator product expansion is given as a series of 5 functions 
and its derivatives and thus is not physical. However, it can be interpreted within the 
shape function formalism. As mentioned before the NRQCD expansion breaks down for 
lepton energies to close to the endpoint. The physical reason for this breakdown comes 
from the lack of phase space for radiating off soft gluons if the lepton energy reaches 
it maximum value. Hence in the endpoint region the non-perturbative subprocess B c — > 
(be) + soft degrees of freedom has to be parametrized by so-called shape functions reflecting 
the phase space dependence rather than just by numbers like the NRQCD matrix elements. 

The formal construction of the shape functions in B c decays proceeds analogously to the 
one in quarkonia decays The phase space of the short distance subprocess is given by 

dPS = dpdp~, t (2^f5\P + l-p t - P u e ) (55) 

where P = M^v and I are the main and the residual cms momenta of the heavy quark- 
antiquark system in the meson rest frame, respectively. If we integrate over d i p Dl and 
expand the lepton energy spectrum in the small residual cms momentum we recognize 
that only the light cone component 1+ = h ■ I is relevant for the soft gluon dynamics. Here 
the light cone vector is given by n = 2(P — p^/M^- Then the lepton energy spectrum 
readsQ 

= ^[n] [dl + f B M(h) S(Mbc/2 -E e + l + /2) (56) 

where the shape functions are defined via 

f Bc [n](l+) = (B c \ [^JC nX ] 8(1+ - iD + ) [xXv] \B C ) . (57) 



As for the matrix elements there is one shape function per each (be) configuration indicated 
by the kernels K, n and fC' n in Q57p. In case of spin singlet configurations vacuum saturation 
is applicable and the covariant derivative iD + reduces to the time component which picks 
up the binding energy M — M&g of the B c . Thus the shape function shifts the endpoint of 
the lepton energy spectrum from the unphysical partonic value M;, 5 /2 to the physical M/2 
expressed by the hadronic B c mass. 

Since the shape function resums purely kinematical effects originating from the phase space 
the short distance coefficients C[n] in eq. (|56|) are not affected at all. This means particularly 
that the coefficient ( pE5|) of the shape function /bJ 1 ^ ] is still helicity suppressed while the 
coefficient ( pH| ) of the spin triplet shape function /b c [ 3 5'| 1 ' ) ] do not vanish for e —>■ 0. 

After changing to our four- dimensional notation and expanding the shape functions in 
the residual momentum iD+jM^ we reobtain the spectrum in terms of delta functions 

2 We ignore the lepton mass mg in this discussion. 
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and their derivatives like eq. (^)- Acting on the spin triplet structure of the dimension- 
6 operator On x i)^S[ ] the zeroth component of the covariant derivative generates the 
following dimension-8 operators (cf. eq. (|47D): 



h<i+) ^ hi-) 

'% ILL "'V 



M, 



be 



h.c. 



Am 2 



l-rf 



3 C (l)l 



r3c(i)i 



(5f 



We re-observe the same structure as we have seen in the pieces of the coefficients ( |52d| - 
f)2f j that remain in the limit e — * 0. Actually these pieces are completely determined by 
the coefficient (@D of £> (lxl) [ 3 S{ 1] ] and the expansion of the corresponding shape function. 
Similarly one can derive the non- vanishing pieces in ( |54b| ) and ( |54cj ) by applying iD±/Mf, 5 
two times on 0( lxl )[ 3 S'j 1 ' ) ]. 

Note that the shape function do not determine the complete coefficients B\ and £> 2 . The 
pieces proportional to e contain terms which cannot be directly derived from the shape 
function. Hence the structure of the spin singlet coefficients in B\ and £>2 is more compli- 
cated. Furthermore we want to stress that the origin of the contributions to Bq which do 
not vanish in the limit of massless leptons is not kinematical. To explain these pieces it 
is not sufficient to introduce shape functions. Rather one has to consider terms that are 
generated dynamically by the non-perturbative physics inside the B c taken into account in 
the bilinear expansion (p^|). 



3.2 Moments of the spectrum 

Another way to understand the unphysical series of delta functions and their derivatives 



is by the means of moments which implies some smearing |22|] 



The two-body structure of the decay at tree level yields a spectrum in terms of an expansion 
in singular functions 

k=0 

where the coefficients of the expansion are given by the moments 

l+e 

M k = [dz^(l + e- z) k . (60) 
./ dz 



A comparison of the result (|39|) with data can thus be performed by sampling the data 
into moments, for which we get a theoretical prediction. The zeroth moment Mq is simply 
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the total rate for which we get 



A + 



(2m) 2 

while the first an the second moments are at least of order 1/m 2 



Mi 



(2m) 



Mo 



r n 



(2m) 



Bo . 



(61) 



(62) 



Hence measuring the second moments of the lepton energy spectra in the semi-leptonic B c 
decay provides direct access on the order of magnitude of the RCM x RCM sector. The 
fact that the values for the matrix elements do not depend on the type of the lepton in the 
final state may help to improve the accuracy of this determination. Likewise measurements 
of the first and zeroth moments may provide some information for the matrix elements of 
the other operators. 



3.3 Total rate in the limit mi — > 



Finally we discuss the limit of vanishing lepton mass. In this limit only subleading terms 
in v survive since the leading term has to vanish due to helicity arguments. On the other 
hand, this limit is certainly a good approximation in the case of the electron spectrum. 

For e — > non-vanishing contribution to the lepton energy spectrum fl39|) contains a piece 
which is proportional to 8(1 — z): 

dT „ f /2 



dz 



(B c \0 {lxl) [ 3 S[ 1) }\B c ) + 



Qm 2 



(6 - 8r] 2 + 3/7 4 ) (O 



15(1)1 



(CxC)l^o ]) 



+ 477(2 - if 



l- v 2 



+ 



12m 2 

(l-?7 
60m 2 



(©(CXR)^]) 



V (o iRxR) M^) + 12 <0 ( r x r)[ 3 p o (1) ]> ] \ 



|5(1 



3 C (l)l 



2\2 



7 <0 ( cxc)flsf >]> - 4 (V^c)^']) 5"(1 - z) 



3 c (l)l 



(63) 



This piece will contribute to the total rate of the semi-leptonic decay B c — > X£p£ even 
if the energy and momentum of the light hadrons in the final state are extremely soft. 
Independent of the exact values of the counting parameters k, A, and v this piece should 
be sizable whereas the purely leptonic B c decay rate ([!]) vanishes for nig — > 0. Consequently 
we conclude that the semi-leptonic rate exceeds the leptonic one by orders of magnitude 
for I = e while it at least reaches it in the muon case (cf. eq. 45). Since we expect the 
light hadrons in the final state to be rather soft the experimental search for the exclusive 
channel B c —> iv^ become very difficult unless the matrix elements in (|63|) are unnaturally 
small. 
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4 Conclusions 



We have investigated the semi-leptonic decay B c — > Xiug where X denotes light hadrons 
below the D° threshold. This decay provides a unique testing ground in quarkonia physics, 
since the two heavy quarks have to annihilate through weak interaction and hence one 
may investigate in this way the light degrees of freedom in a quarkonium. The total semi- 
leptonic rate into light hadrons does not show any helicity suppression at subleading order 
and thus can exceed the purely leptonic one for light leptons, i.e. for electrons and possibly 
also for muons. Thus the measurement of the leptonic B c decay rate could become a quite 
complicated maybe even infeasible experimental task. 

The B c is an intermediate case between the T and the J/if>. While NRQCD is believed to be 
a reasonable approximation at least for the low-lying states of the T, the non-relativistic 
approximation becomes questionable for the J/ip. This point is of relevance for power 
counting, since the usual NRQCD power counting may not be valid for the J/ip- This will 
also affect the contributions related to the centre-of-mass movement of the heavy quark- 
antiquark pair inside a quarkonia-like bound state. For this reason we propose a slightly 
generalized power counting scheme with respect to standard NRQCD velocity rules taking 
care of possible contributions from the cms-motion. 

We have shown that the leptonic energy spectrum provides a testing ground for estimating 
the relevance of such contributions. Measurements of the moments of the spectrum may 
also help to clarify the power counting rules for spin-flip transitions which are a recent 
subject of investigation in charmonia physics. 
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A Bilinears 

To expand the currents T we use the field expansion that can be read off if 
the effective Lagrangian (§) is compared with the full QCD Lagrangian. Additionally 
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considering eq. flTT| ) and the equations of motion ( ]T0| ) this yields 



(64) 



Instead of using the masses of the quark (m+) and the antiquark (m_) we present the ex- 
pansion in terms of the effective mass m and the dimensionless parameter rj parameterizing 
the mass difference: 



2roj.m_ 



m 



and r\ :- 



m_L — rri- 



m + + m- 



m + + m_ 

Then the bilinear expansion up to 0(l/m 2 ) reads: 

QWrgH = &i+>r/£-> 

+ ^ { i&L (hl +) [r, 7,] + hM {r , $ ± } 
+ 17 K (hi +) {r, 7,} ^) + /^ +) [r, <W ^ 



+ (^) {(W) (^(^r/zHJ + ^r^) 2 ^- 



i0£z^ (^ +) -{7^7, + 7,r 7M } hi~ 



I — 

vnr 



(65) 



(66) 



In the case of colour octet case 

• insert T a into the currents hi + ^ . . . hi 

• replace id^_ by the covariant derivative iD^_ 

Finally the corresponding formula for Qi ^ TC is obtained by the replacement rj — > —77. 
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Figure V. Scaling rules for B c matrix elements of operators up to dimension-8. The operator 
power counting is denoted inside the central circle, the Fock state suppression above. 
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